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Abstract. The Spitzer Space Telescope allows for the first time to search systematically for very low luminosity (< 0.1 Lq) 
objects (VeLLOs) associated with dense molecular cores. They may be the first candidate Class sources with sub-stellar 
masses. We describe such a source in the dense molecular core LI 148. VeLLO natal cores show properties that are unusual 
for star-forming cores. The low luminosity and in some cases the lack of prominent outflow could be the result of the small 
gas supply near the VeLLO. 

Key words: stars: formation — ISM: globules — ISM: individual (LI 148) 



©0000 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 

1. Introduction 

The "revealed" phase of early stellar evolution was dis- 
covered and classified in the late 1980's (Lada 1987), but 
it took some years until objects in the "obscured" phase 
of star-formation were found and classified (Andre, Ward- 
Thompson & Barsony 1993). Such Class protostars are so 
deeply embedded in their natal dense cores that they could 
not be detected with contemporary near- and mid-infrared de- 
tectors. They are believed to be in the main accretion phase 
of star formation (Andre, Ward-Thompson & Barsony 2000), 
drive collimated bipolar outflows (Andre et al. 2000) and can 
show maser emission (Furuya et al. 2001). 

While the first Class sources were studied, also the first 
brown dwarfs (BDs) and extrasolar planets were discovered 
(Nakajima et al. 1995). Many BDs are not orbiting a stel- 
lar mass primary (Oppenheimer, KuLkarni & Stauffer 2000). 
There is no proof that these BDs were ejected from form- 
ing protostellar systems. The large BD binary fraction of 
> 15% (Martin et al. 2003, Pinfield et al. 2003) is further- 
more hardly explained by ejection scenarios (Whitworth & 
Goodwin 2005). It thus appears plausible that at least some 
of the single brown dwarfs may have formed by processes 
similar to those by which normal stars form, i.e. the collapse 
of a dense core onto a compact object. Sub-stellar mass com- 
pact objects in the process of formation would appear as low- 
luminosity (iboi ^ 1 Lq) Class sources. 
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It is a challenge to search and study accreting very low 
luminosity objects (VeLLOs) that may turn into sub-stellar 
mass brown dwarfs. While a first such source (IRAM04191) 
was found serendipitously (Andre, Motte & Bacmann 1999), 
a systematic search for such objects has become possible only 
with the Spitzer Space Telescope (SST). Our SST Legacy 
Program "From Molecular Cores to Planet Forming Disks" 
(c2d; Evans et al. 2003), which images about 90 isolated 
dense cores and 5 large molecular clouds (Huard et al., this 
volume), is tailored to find such VeLLOs. The source L1014- 
IRS (Young et al. 2004) was the first low luminosity object 
discovered in our survey. We here present another candidate 
for an accreting sub-stellar object, LI 148-IRS. 

2. The VeLLO Candidate L1148-IRS 
2.1. The compact IR Source 

The dense core LI 148 is part of the Cepheus flare at a dis- 
tance of about 325 pc (Straizys et al. 1992). We have imaged 
the thermal dust emission of LI 148 at 1.2 mm wavelength 
using the Max-Planck Millimeter Bolometer (MAMBO-2) 
aiTay camera at the IRAM 30m-telescope on Pico Veleta 
(Spain). Such wide band continuum observations map the 
projected mass distribution, since the dust emission is nearly 
proportional to the H2 column density and the dust tempera- 
ture (e.g., Motte, Andre & Neri 1998). In our mm-continuum 
map we found two parallel, faint filaments in the south-east 
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2. 1 The compact IR Source 
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Fig. 1. MAMBO dust emission (top) and 24 /im Spitzer map 
of LI 148 (bottom). The bright Spitzer source L1148-IRS is 
marked by a star in the MAMBO map. Offsets are in arcsec- 
onds. 

and north-west of LI 148 (Fig. 0. Assuming dust tempera- 
tures of 10 K we estimate their masses to 1 1 and 5 Mq and 
their peak visual extinctions to 7 and 5 mag, respectively. 

The dust emission peak in the north-western filament is 
associated with a bright Spitzer source at 24 fim and 70 /im 
(Fig.[lli, named LI 148-lRS. From its appearance in the IRAS 
Faint Source Catalog Kun (1998) noted it as a potential low- 
luminosity protostar A spatial, physical association of the 
dust emission peak and LI 148-lRS is suggested by the 6" 
proximity. The probability of a chance alignment within 20" 
of the dust peak (about twice the separation of dust peak and 
LI 148-lRS) and a background galaxy of comparable infrared 
brightness is very low (1:7700 at 24 fim, 1:3400 at 70 /im). 
Because LI 148 is way above the galactic plane (z « 80 pc 
at & = +15°) a chance alignment with a background proto- 
star is also unlikely. Therefore LI 148-lRS is most probably 
an infrared source associated with the LI 148 dense core. 

Using the IRAM 30m-telescope we have unsuccessfully 
searched for a CO outflow from LI 148-lRS, which would 
have clearly identified LI 148-IRS as an accreting protostar. 
This does not exclude the presence of a small-scale flow, 
which in the VeLLO L1014-1RS could only be detected with 
an interferometer (SeclSjl. 
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Fig. 2. The spectral energy distribution of LI 148-lRS. The 
dashed line shows our model SED for the compact emission 
from an embedded protostar (which does not apply to the 
MAMBO data point). 

With the available archival 2MASS and IRAS data and 
our own OMEGA-Cass. (Calar Alto 3.5m-telescope) near-lR 
observations, in Fig.|2lwe plot the observed spectral energy 
distribution (SED) of LI 148-lRS. For A < 100 /im the source 
LI 148-lRS is compact (smaller than a few arcsec) and the 
fluxes listed hold for this compact source. For A > 100 /im 
the source is extended (several arcmin at A = 1 200 fim) and 
the fluxes hold for an aperture of 12'.'9 = 4 200 AU radius, to 
be consistent with previous work (Andre et al. (1999). 

To estimate the bolometric luminosity of the embedded 
star we integrated the flux density over the observed SED us- 
ing a modified version (at A > 200 fim we assume the SED to 
follow a greybody of temperature 10 to 15 K and dust opacity 
K oc A^^ instead of a powerlaw) of the method of Chen et al. 
(1995). The integration yields Lboi = 0.1 to 0.2 (out of 
which about 30 to 70% are due to heating by the interstellar 
radiation field [ISRF]), a submillimeter-to-bolometric lumi- 
nosity ratio L>350/jm/ -^boi > 0.05, and a bolometric temper- 
ature Tboi < 140 K. As L>350Mm/^boi > 0.005 LI 148-lRS 
qualifies as a Class protostar (Andre et al. 2000). 

We tried to model the observed SED with an embedded 
young stellar object (Fig.|3. The continuum radiative trans- 
fer was calculated' for a star of luminosity L and temperature 
T, surrounded by a spherical core with a power-law density 
profile: g cx r^" with a — 3/2 to 2 (Shu 1977) between ra- 
dius Tin and rtrans, a = between rtrans and rout, £» = for 
r < Tin and r > rout- While Tin and rtrans are free parame- 
ters, the outer radius, rout, and density are constrained by the 
extended dust emission seen in the MAMBO maps. The core 
is exposed to the ISRF. We take the differential observing 
technique (e.g., a "background intensity" is subtracted from a 
Spitzer image before a flux measurement) into account when 
modelling the individual flux densities. This produces differ- 
ent model SEDs for the compact (< 2 arcsec at A < 100 /im) 
and the extended emission (> 2 arcsec at A > 100 fim). The 

' The code was kindly provided by E. Kriigel. 



2.2 The Natal Dense Core 
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main result from fitting our models is that independent of our 
choice of values for the free parameters, the stellar luminosity 
and temperature can be constrained to L — 0.05 to 0.15^0 
and T = 2 000 to 4 OOOK. 

In a detailed comparison we notice a discrepancy between 
the predicted and observed flux densities in the 8 /im Spitzer 
band (Fig.|2j, which may be dominated by silicate absorption 
and emission. Furthermore, our model fit implies an inner ra- 
dius, Tin, that is much larger than the dust sublimation radius 
(the dust in our model has temperatures < 100 K, where dust 
sublimates at « 1 500 K), possibly suggesting a more com- 
plicated geometry including the presence of a disk and cavi- 
ties. The presence of significant amounts of hot 100 K) 
dust, in which case we would observe a significant slope in 
the IRAC bands, can be excluded. 

Following Young et al. (2004) we derive an upper limit 
to the protostar mass, M*, using the constraint that the accre- 
tion luminosity, Laccr = GA'hM / is smaller than the to- 
tal luminosity, Lacci < L. Here G is the constant of gravity, 
i?* is the protostar radius, and M the accretion rate. Using 
estimates for i?* and M from Young et al. (2004), we find 
AU < 0.1 Mq. 

2.2. The Natal Dense Core 

Our Effelsberg 100m telescope maps of LI 148 in the 2i- 
lo transition of CCS probe the dense core on large scales 
(Fig. O. The velocity-integrated CCS intensity map reveals 
two filaments associated with the dust emission filaments. 
The CCS centroid velocity map shows a jump in velocity 
about 70" south-west of LI 148-IRS. Since the velocity field 
is smooth elsewhere, the velocity-discontinuity suggests two 
neighbouring, possibly overlapping dense cores, separated in 
velocity by 0.1 km s^^. The linewidth is < 0.2 km s^^ in the 
north-western filament, except where the gradient of the cen- 
troid velocity is large, possibly due to overlapping velocity 
components. 

Our molecular line observations with the IRAM 30m- 
telescope probe the dense core on smaller scales. The C^^O 
2-1 emission, combined with the dust emission implies a CO 
depletion of /d = 8 (Crapsi et al. 2005a); the abundance of 
CO in the gas phase appears to be reduced by a factor 8 rel- 
ative to the "standard" abundance, probably through freeze- 
out onto dust grains during the evolution of the core. Some of 
the observed C^^O emission could come from velocity com- 
ponents unrelated to the dust peak, in which case /d for the 
dense core would be underestimated. 

Our map of the N2H+ 1-0 emission shows that the dense 
gas is closely associated with LI 148-IRS (Fig. |4}, further 
supporting the picture that the Spitzer source is physically 
associated with the LI 148 dense core. 



3. The confirmed Protostars IRAM04191 and 
L1014-IRS 

Besides LI 148-IRS, two other low-luminosity protostars 
have been studied in detail. 
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Fig. 3. CCS 2i-lo integrated intensity, centroid velocity, and 
velocity dispersion. Overlayed contours are dust continuum 
intensity. The position of LI 148-IRS is marked by a star. Off- 
sets are arcseconds. 
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Fig. 4. Dense gas near LI 148-IRS as probed in the N2H+ 
(1-0) transition. The position of LI 148-IRS is marked by a 
star Offsets are in seconds of arc. 
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Table 1. Properties of VeLLOs (top) and their natal dense 
cores (bottom) as observed with single dish telescopes. Dense 
cores are compared with criteria for "evolved" cores from 
Crapsi et al. (2005b); N is column density. A plus indicates 
that a criterion is met, a minus that it is not, a question mark 
that the data is insufficient to decide. Data are from Motte & 
Andre (2001), Belloche et al. (2002), Crapsi et al. (2005b), 
and this work. 



Property 


IRAM04191 


L1014 


L1148 


protostar properties: 








prominent outflow 


+ 








0.45 


0.27 


0.11 


dense core properties: 








U > 10 


? 






7V(N2D+)/iV(N2H+) 


? 


+ 


? 


> 0.1 








iV(N2H+) 


+ 






> 8.5 ■ 10^^ cm"^ 








contr. motions, infall 


+ 




+ 



IRAM04191 was the first known protostar with a very 
low luminosity (L « 0.1 Lq) and therefore mass (Andre et 
al. 1999). It was found through its prominent CO outflow, that 
is associated with a VLA cm radio source and a dust emission 
peak. A near- and far-infrared counterpart was recently iden- 
tified in our Spitzer data (Dunham et al. 2005). 

The VeLLO L1014-IRS was discovered through the as- 
sociation of an unusually bright Spitzer source with a dust 
emission peak (Young et al. 2004). It has a low luminosity 
of 0.1 Lq. The physical association of the Spitzer source 
with the L1014 dense core has been established by the de- 
tection of an infrared outflow cone (Huard et al. 2005; see 
also Huard et al., this volume) and the discovery of a small- 
scale CO outflow with the SMA interferometer (Bourke et al. 
2005). 

4. The Physics of Very Low Luminosity Objects 

Compared with previously discovered Class protostars, 
VeLLOs appear to be suiToundedby small amounts of molec- 
ular gas. Motte & Andre (2001) suggest a comparison based 

on the mass within a projected radius of 4200 AU from the 
protostar, Af 4200 Au^ 

In the Class source sample of Motte & Andre (2001), 
IRAM04191 has the smallest inner envelope mass (and lumi- 
nosity), and the masses forL1014 and LI 148 are even smaller 
than that (Tab. It is suggestive that for VeLLOs the low 
luminosity and lack of outflow activity somehow coiTelate 
with the low envelope mass. However, such conclusion is yet 
based on too small a sample. 

VeLLOs are found in cores that were thought not to be 
able to form stars. Crapsi et al. (2005a) suggested criteria 
for properties of starless cores just before the onset of star- 
formation. According to this, cores must be exceptionally 
dense and chemically evolved, e.g., they show high degrees of 
deuteration and depletion. TablefDshows that the star-forming 



cores L1014 and LI 148 meet only some of Crapsi's criteria. 
Since these were put forward on the basis of single-dish ob- 
servations, they may be insufficient to probe the VeLLO natal 
dense cores on those scales relevant for the star-formation 
process. 

We remark (as Andre et al. 1999 did for IRAM04191 and 
Young et al. 2004 did for L1014) that VeLLOs are embedded 
in dense cores with sufficiently large masses that they could 
in principle grow to stellar masses (> CIA/©), and thus do 
not necessarily remain sub-stellar objects. 

5. Summary 

The Spitzer Space Telescope allows for the first time a sys- 
tematic search for very low luminosity (L « 0.1 Lq) embed- 
ded, compact objects (VeLLOs), accreting Class sources of 
apparently sub-stellar mass. With LI 148-IRS we here present 
a further candidate for such an object. 

The presently known VeLLOs have a very low inner en- 
velope mass. This could relate to their low luminosities and 
- in the case of L1014-IRS and LI 148-IRS - to their low 
outflow activity. Two of the three sufficiently studied VeLLO 
natal cores (L1014-IRS and LI 148-IRS) have properties that 
do not meet the proposed criteria for starless cores before the 
onset of star-formation. They appear neither very dense nor 
chemically evolved when probed with single-dish telescopes. 
VeLLOs seem to form in clouds that were not expected to 
form stars. 
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